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In order to study twist-3 and transverse momentum dependent parton distributions, we use light- 
front time-ordered pQCD at order Qs to calculate various distribution functions for a dressed quark 
target. This study enables us to investigate in detail the existing relations between twist-3 and 
transverse momentum dependent parton distributions. Our calculation shows explicitly that two 
^_^ ■ versions of such relations, considered to be equivalent, occur in the literature which need to be 

^^ ' distinguished. Moreover, we examine sum rules for higher twist distributions. While the Burkhardt- 

^N , Cottingham sum rule for 32 is fulfilled, the corresponding sum rule for /i2 is violated. 

^— ' ' I. INTRODUCTION 

> 

<^ . . . 

[~^ . In view of the increasing accuracy of recent and planned high energy scattering experiments, more and more 

^^ ' attention is paid to the study of parton distributions which are of higher twist and (or) dependent on the trans- 
r^ , verse momenta of the partons. The twist-3 distribution functions are accessible through the measurement of certain 
asymmetries in polarized DIS |l|] and Drell-Yan processes ^J. The transverse momentum dependent (fc^^-dependent) 
structure functions play an important role both in Drell-Yan processes and semi-inclusive DIS ^. In such reactions, 
e.g., the transverse momenta and the transverse spin of the partons can get coupled giving rise to azimuthal asym- 
metries (see e.g. Refs. y-|7|l)' ^hich are very suitable observables for studying the correlations of quarks and gluons 
i"' in hadrons. Often, effects due to higher twist and transverse momenta appear simultaneously like in the recent HER- 
C^. MES measurements of the longitudinal single spin asymmetry in semi- inclusive pion production B . In this work, we 
study these higher twist and /cy-dependent structure functions and their interrelations in the framework of light-front 
hamiltonian QCD. 

As is well-known, twist-3 and /cy-dependent parton distributions are related [H,pl&il4l as a consequence of Lorentz 

invariance. These relations impose important constraints on the distribution functions, which allow one to eliminate 

\^ unknown structure functions in favor of known ones whenever applicable. Consequently, they have been used frequently 

in the literature to facilitate matters, for instance in studying the evolution of /cy-dependent distribution functions 

Our motivation here is to investigate the validity of these Lorentz invariance relations by explict calculation of 
all the involved distribution functions. There exists a very convenient tool based on the light-front hamiltonian 
description of composite systems utilizing many-body wave- functions, which enables us to study these relations in the 
context of perturbative QCD. This tool has already been used successfully in the literature to calculate unpolarized 
and polarized parton distributions ||l^ as well as the transversity distribution ||lj] . The simplicity of this approach 
has also been exploited to make a critical examination of the Wandzura-Wilczek relation ||l5|. Here, we use the same 
approach to calculate the higher twist and /cy-dependent parton distributions perturbatively to order as, and then 
study the Lorentz invariance relations by employing a dressed quark target. We demonstrate that there exist two sets 
of relations, although assumed to be the same, are not identical in reality. More precisely, only one set of relations 
is verified for a dressed quark target, whereas the drawback in the other case can be traced back to the absence of 
quark-gluon-quark correlators, which seem to be crucial ingredients in the Lorentz invariance relations in a gauge 
theory. 

Our calculation also gives us the opportunity to investigate the sum rules for the twist-3 distributions 52 and /12. 
The Burkhardt-Cottingham (BC) sum rule for 172 fl^l is satisfied for the dressed quark target, but the corresponding 
sum rule for /12 |17] , |l§[| turns out to be violated. To our best knowledge, the violation of the sum rule for /12 in a 
perturbative treatment is a new observation. 
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The paper is organized as follows: In Sec. y, we outline the definition of the parton distributions relevant for our 
discussion, and give a detailed account of their relations due to Lorentz invariance. In Sec. Ill, first the operators for 
twist-3 distributions are expressed in terms of dynamical fields, and some relevant points regarding the dressed quark 
target are discussed. Then we present our results for the different parton distributions and a detailed investigation of 
the two sets of Lorentz invariance relations. In Sec. IV, we study the sum rules for 172 and /i2 and conclude in Sec.M. 
Some conventions are summarized in an Appendix. 



II. PARTON DISTRIBUTIONS: DEFINITIONS AND THEIR INTERRELATIONS 

In this section we recall the definitions of various parton distributions that already exist in the literature and 
introduce the Lorentz invariance relations among them. We restrict the discussion below to twist-3 structure functions, 
while in the case of fc^-dependent functions we limit ourselves to twist-2 level which is sufficient for our purpose. For 
a complete discussion one should go back to the original references mentioned below. 

To begin with, we specify the correlator ^{x) of two quark fields on the light-front[| in terms of which all the 
structure functions are defined. 
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with k^ =^ xP^ . The target state is characterized by its four- momentum P and the covariant spin vector S (P^ — 
M^, S^ = —1, P ■ S = 0). Note that <I>(x) can easily be made gauge invariant by putting an appropriate gauge link 
between the quark fields. However, since the non-localilty in the operator is only in the longitudinal direction ^~ and 
we shall be working in the light-front gauge (yl+ =0), we can always get rid off the gauge link in Eq. p). 

Now, the parton distributions appear in a general decomposition of the correlator $(a;) where one finds three 
functions at twist-2 and three functions at twist-3 level 
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with n_|_ and n_ being two light like vectors satisfying n+ • n_ = 1. The helicity of the target state is given by A, while 
St = (0, 0, St) represents the transverse spin of the target. Sometimes different notations for twist-2 distributions 
are used in the literature (/i(a;) = q{x), gi{x) = Aq{x), hi{x) = AtQ^x) = Sq{x)). The twist-3 part contains the 
well-known transversely polarized structure function gx, and two chiral-odd distributions e and h^. Note that in 
Eq. (0) we have not considered the so called T-odd parton distributions. 

The structure functions in Eq. ^) are projected out by performing traces of ^{x) with suitable Dirac matrices. 
Using the abbreviation ^t'"' = Tr (<I'r)/2, we give the explicit expressions for those structure functions that are relevant 
for our discussion here. 
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where, like in Eq. (|^), all the correlators are understood to be on the light-front, i.e. C+ = ^t = 0- 
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^Our definition of light-front components of a generic 4- vector as well as further conventions are summarized in the Appendix. 



In a similar way, /cT-dependent parton distributions are defined starting from the following correlation function 
where the non-locality in its operator structure is not only in ^~ but in ^^ as well, 



$,, (x, fc-) = ^ / ^90^ '^''^ <^' ^1 ^^(^^ ^^^^) 1^' ^> 



(7) 



Here we have assumed that in A'^ = gauge together with antisymmetric boundary conditions for the transverse 
gluon field, the gauge link can still be omitted as argued, e.g., in Ref. [g|. In the general decomposition of this 
correlator one naturally finds more distribution functions due to the presence of an extra vector k!^ [g , 
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Here we have shown only the twist-2 part, which is sufhcient for our purpose and have omitted the T-odd functions 
as before. 

Like in the previous case, one projects out the structure functions in Eq. (|^) by performing traces of ^{x, fe^) with 
suitable Dirac matrices. Two of the projections necessary for our discussion are, 

-i= ^[''^^^l {x, k^) = XgiLix, kl) + ^^^ gM^, kl) , (9) 

-^ $[-'^^=1 {x, k^) - SI hM^, kl) + ^ (^A hi^{x, kl) + ^^ hi^ix, kl)^ . (10) 

Note that depending on the target polarization the same projection of <I'(x, fc^) allows one to calculate different 
structure functions. For example, from Eq. (0) we get 5il(x, fey) or giTiXjk^) for the target being polarized in 
longitudinal or transverse direction, respectively. Therefore, Eqs. (uTO) give us a well-defined way to calculate 



giT{x, k^) and hj^]^{x, k^) which are necessary for the subsequent discussion. In what follows we need the /c^-moments 
of these two functions, which are defined as p| 

g^)[x)^ jd?kT^^g,T{x,kl), hf^'\x)^ld'kT^hi^ix,kl). (11) 

Now that we have given all the definitions of relevant structure functions, we are in a position to discuss the existing 
relations among them. These are usually of two kinds - one follows from the QCD equations of motion and the other 
comes as a consequence of Lorentz invariance. Here we are mainly interested in the latter ones which, according to 
Refs. [Hq], read as 

gT{x)^giix) + ^g['^{x), (12) 

hL{x)^h,ix)-^hfl^'\x). (13) 

These relations have been derived from the general Lorentz covariant decomposition of the correlation function $ of 
two quark fields before it is constrained on the light-cone and, hence, they are quite naturally referred to as Lorentz 
invariance relations. On the other hand, a similar relation for gr also attributed to Lorentz-invariance has already 
been proposed in Ref. ^ and extended for Hl in Ref. ||l^. A detailed account on these relations can be found 
in Ref. [|ll|. Using the notations of Ref. ||ll|], with necessary modifications for the conventions that we follow, the 
corresponding relations are given by 

f \ / ^ d -. . /■ , , D(x,x') + D(x' ,x) ,_, ,, 

gT[x) ^ gi{x) + —K{x) + dx' ^, _ ^ ^ (14) 

, , , , , , I d ~, , 1 /■ , ,D{x,x') + bix',x) , , 

fti(x) = h,{x) + --i^(x) + - / dx' j,_^^ . (15) 



Note that here quark-gluon-quark light-front correlators are involved, which depend on two momentum fractions 

k+ 



denoted as a; = -mr, x' — 4q- and the new correlation functions appearing in Eq. (n4h are given by 
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while the ones in Eq. ( |l5| ) are 
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with D{x,x') = i[i:»i(x,x') + D2ix',x)] and D{x,x') = i[i:»i(x,a:') + 1)2(2:', a;)] ■ In principle, Eqs. (|T|,|l|) and 
Eqs. ( [l^jl5| ) (henceforth, referred to as Set-A and Set-B respectively) should contain the same information and, in 
fact, assumed to be identical. But there exists hardly any proof of that. 

These relations are quite remarkable, in particular, since they involve at the same time functions describing longitu- 
dinally and transversely polarized targets and therefore will provide us with a consistency check while comparing data 
for the measured structure functions from different experiments. Moreover, they can be quite useful to predict the 
evolution of one of the structure functions once the evolutions of others are known, as has been done in Refs. P,ttO,n2|. 
Keeping their importance in mind, it is worthwhile to delve more into these relations. We do this in the next section 
by checking them through explicit calculations for a dressed quark target in the framework of light-front time-ordered 
pQCD. 

III. PARTON DISTRIBUTIONS: RESULTS AND DISCUSSION 

Before presenting our results, it is useful to disentangle the twist-3 parton distributions into simpler structures which 
manifest the different aspects of the QCD dynamics contained in them. To achieve this, we re-express the structure 
functions gr and h^ in terms of dynamical fields, the so-called good fields, like in the twist-2 case which right from 
the beginning contains only the good fields. That is, we eliminate the constrained field '0_ via the constraint equation 



(A2) in terms of ■0-1- and A\, which are the only dynamical fields in the hamiltonian formulation of light-front QCD. 



Thus, gx and h^ defined in the previous section become 

= g^{x,Q^)+g^T-{x,Q^)+g?r{x,Q^), (22) 

= hl[:{x,Q^) + hl-{x,Q^)+hi{x,Q^), (23) 

where we have introduced the operators 
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Here ruq is the quark mass and At = X^a-^^^r the transverse gauge field, while 1/9+ is defined in the sense of 
the principal value prescription as given in Eq. (A3). The above light-front expressions make the physical picture 
of twist-3 structure functions clear. It explicitly shows the contributions associated with the quark mass, quark 
transverse momentum and quark-gluon coupling operators. Although one naively expects that the contributions 
depending explicitly on the quark mass are suppressed, it turns out that each of them is equally important to extract 
the information contained in twist-3 structure functions. Notice that Eqs. ( |2^ , |23| ) correspond to what in the literature 
is often referred to as the relations among various light-front correlators coming from the QCD equations of motion 
(see e.g. Ref. |l). 

Having presented the relevant issues as far as the operator structures involved in the parton distributions are con- 
cerned, some comments regarding our calculation and the target state are in order. The calculation is straightforward 
and we shall avoid giving unnecessary details except mentioning the following points. (For details we refer the reader 
to Ref. |l|.) 

Firstly, all the required structure functions are calculated for a dressed quark target given by the following Fock- 
space expansion truncated at the lowest non-trivial order, 

\k, A) = M\b\{km + Y^ f ^fc^^^fc^i dkU'kT2 ^2(^ ^3(fc _ fc^ _ fc^) 
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where &^(fc) and aj^(fc) are the creation operators of quarks and gluons on the light-front which obey the usual 



commutation relations (see Eqs. (A8;A9)). The most important ingredient in the above dressed quark state is the two 



particle boost-invariant wave-function which can be calculated using light-front time-ordered pQCD and is given by 



*U(-,-.) = -^|^^^xl|2^ + -(-.-.)<— .<^|XA.-A., (27) 

with X and Kt being the relative momenta of the quark. Note that the mg-dependence in the above wave function has 
its origin in the helicity flip part of the light-front QCD Hamiltonian. This is an essential term in investigating the 
dynamics of transversely polarized targets and, hence, is also very important as far as our calculation is concerned. 
The constant Af appearing in Eq. ( |26[ ) is determined by the normalization condition 

{k', A'|fc, A) = 2(27r)^ k+ 5{k+ - k'+) S^ik^ - K) (5a,a' , (28) 

and to the order as given by |13t| 
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Here a hadronic scale fi has been introduced such that k/^ >> /i^ >> (m,)^, which can be considered as the 
factorization scale separating the "hard" and "soft" dynamics of QCD. This scale fi also serves as the lower cutoff of 
the involved transverse mometum integration, whereas Q^ is the upper cutoff. 

Secondly, in our calculation we also need a transversely polarized target, for example, in the case of gx- This is 
obtained by a superposition of two different helicity states. Thus, the one polarized in x-direction can be expressed 

by 



|fc,^i = ±i) = -i=(|fc,T)±|fc,i) 



(30) 



Lastly, the quark mass renormalization enters in the calculation at a^-order and we use the following expression 
for the rcnormalized quark mass m^ in terms of its bare mass rUg [^ , 



mf = m,{l+'-^Cjln^ 
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We now present the results of our calculation, i.e., all the relevant structure functions for the dressed quark target 
in Eq. (26) up to order as- We first give the twist-3 structure functions gx and hj^. It turns out that all the three 
terms in Eq. (E3) and Eq. ( p3| ) have nonzero contribution to the corresponding twist-3 structure functions and for 
clarity we provide them separately. For gxi^) we obtain 
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A similar calculation for h^ gives 
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Note that the above results represent purely the pQCD dynamics to the a^-order. It should be noted that all the 
individual contributions in gr as well as h^ in the perturbative calculation are of the same order (namely, proportional 
to niq/AI), which means that the mass dependent terms g™ and /i™ are not suppressed contrary to the common belief. 
As mentioned above, nig is the bare quark mass and up to order as it is given by Eq. ( pl| ) in terms of the rcnormalized 
quark mass m^. On the other hand, M is the rcnormalized target mass and, therefore, in our case it is identical to 



m^ itself, M 



Taking this into account, we finally get gr and h^ as follows. 
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where we have used the well-known 'plus'-prescription. Eq. (|3|) reproduces^ the result already obtained in Ref. |l^] 
for a dressed quark target. Also a covariant one- loop calculation with a quark target yields exactly the same expression 
for gT pG|. Note that the result for ft^, obtained for the first time here, does not contain any singularity at cc = 0, 
even though the individual pieces in Eqs. ( p5[]37| ) do. 

To investigate the validity of the relations in Set-A, we need to calculate the structure functions on the RHS of them 
which involves explicitly calculating gi, gYp^ h\ and h^)^ for the dressed quark target. Carrying out the evaluation 
of g\fl and h-^)^ as given in Eq. ([ll|), we get to the a^-order 



■^Our result differs from that obtained in Ref. |I5| by a factor of i which appears in the definition of qt that we are using 
This is not relevant for our purpose as long as we use one consistent set of definitions for all the parton distributions. 
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We point out that in the free theory [as = 0) these two /cT-dependent functions vanish since our target carries no net 
transverse momentum. The same is true for the two functions given earUer in Eqs. ( |33| , |36| ) . 
The results for gi and hi to a^-order already exist in the hterature |l3,M and are given by 
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Note that the ag-terms for gi and hi contain the evolution kernels of the corresponding structure functions. Having 
the explicit results for all the necessary structure functions appearing in the Lorentz invariance relations as given in 
Set-A, we can now compare the LHS and RHS of these relations. By doing so, one readily finds that the relations in 
Set-A are not satisfied for a dressed quark target. Therefore, the natural conclusion is either that Lorentz invariance 
is violated in perturbation theory or that the relations in Set-A do not reflect the complete picture. 

It is easy to see that these relations in Set-A in fact do not reflect the complete picture. To make it evident, we 
turn our attention now to the relations in Set-B and first calculate K and K for the dressed quark target to the same 
order. It turns out that 
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This immediately leads us to the conclusion that the relations presented in Set-A are actually different from that in 
Set-B unless the contributions coming from the D(a;,a;')'s and Z3(a;,a;')'s are identically zero, which is unlikely in a 
general scenario. Carrying out the explicit calculation of these terms, we get 
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We point out that no singularity in (x' - 
the apparent pole there gets c ancele d. 
Putting the results in Eqs. 



x) shows up in Eqs. (|1^,|15|). From the above results one easily observes that 
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3|) back in the RHS of Eqs. (^_4||1^), we obtain gx and /i^ as given in Eqs. 
which verifies the relations in Set-B. Moreover, we see that the discrepancy we found earlier in Set-A is exactly 
compensated by taking these D{x,x')'s and D{x,x')^s properly into account. In other words, from this excercise it 
turns out that the information contained in these quark-gluon-quark correlators is missing in the relations given in 
Set-A, thereby making them incomplete. 

Therefore, we finally conclude that the relations in Set-A and those in Set-B are not identical - while the first ones 
are violated, the second ones are fulfilled for a dressed quark target up to order as- Barring Eqs. (p^ 34 ) , (38) and 
(E3,E3|), all the results presented in this section are obtained for the first time here in the context of light-front QCD. 
It should be noted that in the free theory both sets of relations are satisfied, which is easily verified by setting all the 
terms proportional to a^ in the above expressions for the structure functions to zero. Since only the quark-gluon-quark 
correlators seem to be missing in Set-A, we believe that it is valid and useful in models where no gauge fields are 
involved. For instance, we have checked explicitly taking the results for the parton distributions as obtained in the 
spectator model [E2i that the relations in Set-A can be verified. However, in the context of a gauge theory like QCD, 
one should be careful and always be using the relations in Set-B. 

IV. SUM RULES 

Our calculation here provides a direct way to investigate the existing sum rules for twist-3 parton distributions in 
the case of a dressed quark target. Defining the structure functions g2 = gr — 3i and /12 = 2 {h^ — hi), the following 
sum rules have been proposed in the literature pq-p^: 
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From the results presented in the previous section (see Eqs. ( P8[|39| ) and (42|43|)) we can immediately write down the 
expressions for 52 and /12 to the order as, 
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The BC sum rule for 52 follows readily from Eq. (48) as has already been shown in Ref. |l5l using the same method. 
On the other hand, Eq. (49) gives 



dx h2(x) = ^- Cf In — ; 



(50) 



which shows that the /12 sum rule is violated in perturbation theory. Incidentally, the second moment of /i2 turns out 
to be zero, /„ dx xh2{x) — 0, although the significance of this result is not clear. 

The violation of the /i2 sum rule in the context of perturbation theory is a new observation. Such a result is 
unexpected, bearing in mind that the ft,2 sum rule has been derived on the same footing of rotational invariance like 
the BC sum rule [^. Wc point out that our observation is different from the findings outlined in Ref. ||2^, where 
a possible violation of this sum rule for the experimentally measured structure function has been discussed. There, 
the origin of such a violation was attributed to quark zero modes giving rise to a (5-function singularity in the parton 
distribution at x = 0. Since this kinematical point is usually inaccessible, a significant deviation from the sum rule 
could occur in the experiment. In contrast, we find the violation already at the level of the parton distribution /12 
calculated to the order a^. Moreover, in our explicit calculation the final result for ft,2 is not inflicted by quark zero 
modes. 



V. CONCLUSIONS 



In this work we have calculated higher twist and fc^-dependent parton distributions using the light-front hamiltonian 
description of composite systems in terms of multi-parton wave- functions. Employing a dressed quark target we have 
evaluated them to the order as in light-front time-ordered pQCD. While we have reproduced the results for gi, hi 
and gT, all the other results presented in Sec. |l| are new. 

These calculations, in particular, have given us the opportunity to study the so-called Lorentz invariance relations 
existing among twist-2, twist-3 and fcr-dependent structure functions. We show explicitly that two distinct sets of 
such relations exist in the literature. While one set is satisfied (Eqs. (n4,n5h) for the dressed quark target, the other 
one (Eqs. (02,0)) is not. It turns out that quark- gluon-quark correlators arc important for the Lorentz invariance 
relations, where these pieces are what exactly is missing in Eqs. (0,13|). The implication of our findings on the 
existing literature is yet to be explored. 

Moreover, we have studied the sum rules for the structure functions 52 and /i2. The BC sum rule for (72 is fulfilled, 
whereas the corresponding sum rule for ft,2 is violated at the order as in perturbation theory. Since both sum rules 
have been derived on the same basis of rotational invariance, the violation of the /i2 sum rule is surprising and requires 
further investigation. 
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APPENDIX A: 



In this Appendix wc summarize our conventions. First, we specify the plus and minus hghtcone components 
of a generic 4- vector a^^ according to a^ = a^ :k a?' ^ and the inner product of two 4- vectors is given by a • 6 = 
ia+6^ + \a^b^ a^ ■ b^ . For the 7 matrices we use the hght-front representation p5[ 



7 



-i 
i 



7 = 



-i (Ji 
idi 



1' = 



i 
i 



7^ = 



(TZ 

-fT3 



(Al) 



where di = a2 and (T2 — — cri. In the usual way we define the dynamical field ■0+ = A+V' and the constrained field 
'0_ = A"?/!, which follows the constraint equation 



7 

where Z?^ = d^ — igs^'^ is the covariant derivative. The operator 1/9+ is defined as 

1 1 f°° 



(A2) 



(A3) 



with e{x) being the sign-function. In the representation (Al), the projection operators A^ = 7+7^/4 take the simple 
form 



For the fermion fields we use the two-component notation p2 




1 



^+ = I S ) . ^- 



j^i^T ■ {idr + 9sA-r) + imq\ri J ' 



where the Fourier expansions of the dynamical fields 'ri(x) and Aj,(x) are given by 






2(27r)3Vl+ 



bxik)e 



-ik-x , jf 



dlxik)e 



ik-x 



\ax{k)eT. 



\e 



-ik-x j^ t 



4(fc)4,A 



Ak-x 



(A4) 
(A5) 

(A6) 

(A7) 



2(27r)3fc+ 

Here the creation and annihilation operators for quarks (gluons) obey the anticommutation (commutation) relations 

[bx{k), b{,{k')] = {d;,(fc), 4,(fc')} - 2(27r)3 k+ 6{k+ - k'+) 5{K - k'r) Sx,x' , (A8) 

ax{k), al, (fc')l = 2(27r)3 k+ S{k+ - k'+) S{kr - K) -^a.a' • (A9) 
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